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Calculating the Isgur-Wise Function on the Lattice 
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We calculate the Isgur-Wise function by measuring the heavy-heavy meson transition matrix element on the 
lattice. The standard Wilson action is used for both the heavy and light quarks. Our first numerical results are 
presented. 



1. Introduction 

Because of the new flavor and spin symme- 
tries in the heavy quark effective field theory 
(HQEFT), the heavy-heavy meson decay matrix 
element can be simplified and different decay pro- 
cesses can be related to each other. For example, 
in the case of B — > D decay the number of un- 
known form factors can be reduced from two to 
one, and we have Q 

< D v i\cry„b\B v >= ^/m B m D C c i>£,{v' -v)(v+v') v , (1) 

where v v , v' v are the four- velocity and tub , tud are 
the B and D meson mass, respectively. The con- 
stant C c b comes from integrating the full QCD 
contribution from the heavy quark mass scale 
down to a renormalization scale \x <C 
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where a(v ■ v') is a slowly varying function of v ■ v' 
and vanishes at v = v' The Isgur-Wise func- 
tion £(v'-v) represents the interactions of the light 
degrees of freedom in the heavy meson system and 
can thus be calculated only by nonperturbative 
methods. 

On the lattice the heavy meson system can be 
studied in two different approaches. One is to 
keep the heavy quark dynamical by using the 
standard Wilson action. This may require ex- 
trapolation to the physical heavy meson mass of 
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interest. The alternative is to integrate out the 
heavy quark first and derive an effective action in- 
cluding only the light degrees of freedom and then 
perform numerical simulation using this effective 
action. Here we stay with the first approach. 

Using the flavor symmetry of HQEFT the 
Isgur-Wise function relevant to the B — > D decay 
of Eq. ([!]) can be obtained also from the D — > D 
elastic scattering matrix element |l| 

< D v i\crf v c\D v >= m D C cc (n)^{v-v'){v J rv') u , (3) 

where 



c cc (p) 



a s {m c ) 



a(v-v') 



(4) 



This of course requires that the D meson be suf- 
ficiently heavy for the onset of the heavy quark 
limit (HQL). Conventionally a B — > D (here we 
use B and D as generic names for heavy pseu- 
doscalar mesons, they do not necessarily repre- 
sent the physical B and D mesons) transition ma- 
trix can be parametrized as 

< D P ,\V V \B P >= f + (q 2 )(p'+P),+f-(<l 2 )(p-p')v . (5) 

where q 2 = (p' — p) 2 is the momentum transfer 
between the initial and final states and V v is a 
vector current. It is easy to show that in the 
elastic scattering case one has f-{q 2 ) = 0. Using 
the relation v v = p v /m, Eq. (|^) becomes 



< D V ,\V V \D V >= m D f+(q 2 )(v' + v\ 



(6) 
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Comparing this with Eq. (J3J) one finds the simple 
relation between / + and £ 



(7) 



The lattice calculation method for / + has been 
well established f|,f§, and thus the result can be 
easily used to obtain £. 

2. Considerations in the Lattice Calcula- 
tion 

At j3 = 6.0 the inverse lattice spacing a -1 w 
2.0GeV. The HQL becomes valid when the heavy 
quark Q in a heavy meson has a mass mg 
Aqcd ~ Q-2GeV. The heaviest mass we can take 
is order of one or less in lattice units, beyond 
which one would expect large lattice-spacing ar- 
tifacts. At j3 = 6.0 this corresponds to a physi- 
cal mass in the range of D meson. This can be 
obtained by setting the hopping parameter for 
the heavy quark Q to || kq = 0.118. For the 
light quark we take the hopping parameter n q — 
0.152 — 0.155 and extrapolate to n q ^ cr — 0.157 ||. 

How far can v ■ v' change on the lattice? In our 
calculations we always have either the initial or 
the final particle at rest. Thus 



E D /m D if v' = (0,0,0,1) 
E' D /m D if v= (0,0,0,1) 



(8) 
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where Ed = y m \ 

For a spatial lattice size L — 24, we have in- 
jected momenta 



_ 2tt . , 2ir , 

p= —(1,0,0) ,—(1,1,0) 



(9) 



Since mo w 1.0 in lattice unit, one gets v ■ v' « 
1.034 and 1.066 respectively. To get larger values 
for v ■ v' one needs to inject larger lattice mo- 
menta which would in turn introduce large sta- 
tistical noise in the matrix element calculations. 
Thus in practice v ■ v' can not be much more than 
- 1.1. 

However, this apparent restriction in the range 
of v-v' in the lattice calculations has little physical 
consequence as the validity of HQL requires that 
the momentum transfer between the initial and 
final light degrees of freedom, ~ ^qcd( v ' v ' — 1)> 



be <C Aq C£) [|IJ. This in turn means that v ■ v' 
should be close to one. 

Removing the lattice artifacts. One of the ma- 
jor concerns in this calculation is the size of the 
lattice artifacts. Since the heavy meson mass is 
near one in lattice units, the lattice artifacts could 
be significant. Ultimately the lattice artifacts 
can be brought under control either by compar- 
ing data at different (3 values or using the lattice 
improved actions. However, for simulations at a 
given P value there are several ways to check the 
size of the lattice artifacts. 

Eq. (|^) and consequently Eq. (|6|) require 
Lorentz invariance to hold. In Euclidean space 
the Lorentz transformation becomes a four di- 
mentional Euclidean rotation. The lattice the- 
ory, however, does not have exact Euclidean rota- 
tional invariance for finite lattice spacing a. Thus 
/_ is not exactly zero. The amplitude of /_ (or 
/_//_!_) gives a measure for the violation of the 
Euclidean invariance on lattice. 

We can also estimate the size of the lattice ar- 
tifacts by checking the simulation results against 
known continuum matrix element values at some 
special points. For example, when both the ini- 
tial and final D mesons are at rest, v = v' = 
(0, 0, 0, 1), the continuum matrix element of C74C 
is known because of the quark flavor current con- 
servation m 



< D\c^4c\D >= 2m D 



(10) 



At this so-called "recoil-point" we have £(1) = 1. 
Both Eqs. (jl^) and (j|) will have O(a) corrections 
on lattice. These lattice artifacts may come from 
different origins. Part of the 0(a) effect can be 
approximately included by using a normalization 
factor 

< ip(x)tp(0) > CO nt= 2nu e m < ip{x)$(0) >i att ,(11) 
where [01 



e ™ = l + ifJ--J- 

uq \2k 2k c 



(12) 



with uq the "tadpole improvement" factor. 

Another correction comes from the use of (non- 
conserved) local vector current V v — cj v c ||. 
This effect can be corrected by introducing a 
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Table 1 

The Isgur-Wise function. The heavy quark hopping parameter kq and the lattice sizes are shown in the 

table. 

.118, 16 3 x 39 .118, 24 3 x 39 .135, 24 3 x 39 

1.0567(8) 1.1103(15) 1.0259(5) 1.0512(10) 1.0543(12) 1.1059(23) 
£ 1.00(4) 0.99(11) 0.974(20) 0.940(50) 0.952(20) 0.893(40) 



rescaling factor Zy c in the vector current. In per- 



turbation theory Zy c is calculated to be 



1 - 27.5- 



16tt 2 



(13) 



In general the size of lattice artifacts will be 
momentum dependent so the 0(a) correction will 
be different for Eq. (10) and Eq. (^). However, 
since we are using only small momentum injec- 
tions we may assume that the leading 0(a) cor- 
rection does not depend on p, p' and the Lorentz 
index v. We will use Eq. ( |j"o| ) as the normaliza- 
tion condition for matrix element < D v > | V v \ D v > . 
This way both exp(m) and Zy c factors are taken 
out automatically (together with any other mul- 
tiplicative O(a) factors). 

3. Numerical results 

We use the standard Wilson action for quarks 
in the quenched limit. Both heavy and light 
quarks are treated as dynamical. We use data at 
(3 = 6.0 on 16 3 x 39 and 24 3 x 39 lattices. There 
are 19 configurations on the 16 3 x 39 lattice and 8 
configurations on the 24 3 x 39 lattice. The tech- 
niques for measuring the two-point function and 
the three-point matrix elements are standard || . 
The fittings are done for time slices from t=10 
to t=15. To reduce the statistical error, we have 
used symmetry properties of the Green functions 
and averaged over ±t and ±p directions. 

Comparing the measured / value at q 2 = 
(/o(0) = f+(q 2 = 0)) to the known continuum 
value /o(0) = 1, we observed that the lattice ar- 
tifacts are typically 20% - 40% at k q = 0.118 
and less than 10% at kq = 0.135. We also mea- 
sured the ratio /-//+. We find the violation of 
Euclidean invariance is typically 5% — 15% for 
kq = 0.118 and 3% - 10% for k q = 0.135. Note 
that the reduction of the lattice artifacts when 
kq is changed from 0.118 to 0.135 agrees with 



our intuitive expectations. One may try to use 
the factors Zy c and e m to remove part of the 
0(a) effects. At (3 = 6.0 we get from perturba- 
tive calculation Eq. @ that Z l {/ C w 0.7 (using 
the shifted effective gauge coupling g 2 w 1.7 as 
suggested in ref ^|) . The factor e m is about 2 for 
kq = .118 and 1.5 for kq = .135. Including both 
these factors the corrected /o(0) becomes 1 within 
errors. This is in agreement with other observa- 
tions II that Zy c and e m factors seem to account 
for the largest part of the lattice artifacts. 

According to the discussions in Section 2, we 
define £(v ■ v') = /+//o(0) and list the results 
in Table 1. We emphasize that this definition 
removes all momentum-independent O(a) effects 
simultaneously, including both e m and Zy c fac- 
tors. Note that in Eq. (^) there is a factor C cc 
in the connection between / + and £. This factor 
comes from integrating out the QCD effects from 
the heavy quark scale down to a light scale /i. For 
the lattice calculation, however, fi is taken to be 
0(1/ a). Thus in our case /i ~ mjj ~ m c . Also 
vv' is very close to one. So for practical purposes 
we can set C cc = 1 according to Eq. ([!]). 

We plot our results for the Isgur-Wise function 
in Fig. 1. For comparison we also plotted the the- 
oretical bounds on the Isgur-Wise function. The 
top and bottom curves are the upper and lower 
bounds derived in ref. |J. They are obtained 
using the dispersion relation for the two-point 
functions, with the requirements of unitarity and 
causality and with some assumptions on the ana- 
lytic properties of the form factors. The curve in 
the middle is an upper bound on the Isgur-Wise 
function derived from current-algebraic sum rules 
]Tc| ]. This is a tighter upper bound. Our data ob- 
tained on the 24 3 x 39 lattice appear to be, within 
errors, inside of the upper bound of Bjorken |l0[ 
and the lower bound of de Rafael and Taronpf. 
The data from the 16 3 x 39 lattice is consistent 
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Table 2 

Comparison of the lattice calculation for the slope at v ■ v' — 1 with various model calculations. 
Lattice BUB B B§ 

1.0(8) 1.6(4) 1.4(6) 1.05(20) 0.65(15) 0.25 < p z < L42~ 



with the Bjorken upper bound within rather large 
errors. 



Figure 1. The Isgur-Wise function is plotted 
against v ■ v' . The open circles represent data 
on the 16 3 x 39 lattice. The solid circles and 
squares represent data on the 24 3 x 39 lattice. 
The heavy quark hopping parameter is set at 
hq = 0.118 (open and solid circles) and at k = 
0.135 (squares) while the light quark hopping 
parameters are extrapolated to the chiral limit 
k = k c ~ 0.157. 



Close to v ■ v' = 1 the Isgur-Wise function can 
be parametrized as 

£(v ■ v') = 1 - p 2 {v ■ v' - 1) + 0((v ■ v' - l) 2 ) .(14) 

If we calculate the slope using the data point clos- 
est to the v ■ v' = 1 axis, we get p 2 — 1.0(8). In 
Table 2 we list our result along with p 2 values 
estimated by other authors. 

Note that the lattice meson mass we used is 
in the range of physical D meson. Thus the 
0(l/rriQ) correction may be quite significant. For 
a reliable calculation of the Isgur-Wise function 
we need to repeat this calculation for several dif- 
ferent masses and then extrapolate to the infinite 



mass limit. For a check on the residual O(a) ef- 
fects we plan to repeat our calculation at f3 = 6.3. 
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